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Abstract 

We consider the following geometric optics problem: Construct a 
system of two reflectors which transforms a spherical wavefront gener- 
ated by a point source into a beam of parallel rays. This beam has 
a prescribed intensity distribution. We give a rigorous analysis of this 
problem. The reflectors we construct are (parts of) the boundaries of 
convex sets. We prove existence of solutions for a large class of input 
data and give a uniqueness result. To the author's knowledge, this is 
the first time that a rigorous mathematical analysis of this problem is 
given. The approach is based on optimal transportation theory. It yields 
a practical algorithm for finding the reflectors. Namely, the problem is 
equivalent to a constrained linear optimization problem. 



1 Introduction 

We consider the following beam shaping problem from geometric optics, sketched 
in Figure 1: Suppose we are given a spherical wavefront with a given intensity 
distribution emitted from a point source. We would like to transform this 
input beam into an output beam of parallel light rays with a desired inten- 
sity distribution. This transformation is to be achieved with a system of two 
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Figure 1: Geometry of the reflector problem. Note our convention regarding 
spatial coordinates, illustrated by the coordinate system in the lower left hand 

corner: The output beam propagates in the direction of the negative 2;— axis, 
and points in the plane perpendicular to the 2;— axis are denoted by the vector 
X e M^. Thus a generic point in three dimensional space is denoted by {x, z). 
For more details see the accompanying text in Section 2. 

reflectors^. This paper deals with the mathematical problem of flnding these 
reflecting surfaces for given input and output apertures and input and output 
intensities^. The geometric optics approximation is assumed. 

These types of problems are of practical interest and there exists an exten- 
sive engineering literature on them; see Section 5 in V. Oliker's stimulating 
survey paper [6] and the many sources cited there. As explained there, the ap- 
proaches described in the engineering literature are usually only applicable to 
specific data and are usually only justified by specific numerical examples. In 
[6], the problem was reformulated as equation of Monge- Ampere type for the 
polar radius of the first reflector. A rigorous mathematical analysis, showing 
the existence of solutions, was however lacking as stated in [6] . 

In the present paper, we provide this rigorous mathematical analysis, using 
a novel approach to the problem. (To the best of our knowledge, this is the 
first time such a rigorous analysis has been provided.) We show the existence 
of solutions for given input and output intensities and provide a uniqueness 
result for the ray tracing map associated with the reflector system. The two 
reflectors we obtain are always concave; they are contained in the boundary 
of certain convex sets. The approach is based on advances in the theory of 
Monge-Ampcre equations and optimal transportation that have been made 
in about the last 15 years [1, 2, 3, 10]. Similar approaches have fruitfully 
previously been applied to other beam shaping problems by X- J Wang in [12] , 
and independently by V. Oliker and the author in [5, 4]. 

We now describe the approach and results in more detail. For this, denote 
the input aperture by and the output aperture by T. Thus D is & subset of 

-'^A treatment with three reflectors is possible, sec [8]. For reasons of space limitations, a 
two reflector system is sometimes more desirable from an engineering point of view. 

^From a practical point of view, an additional requirement is that no blockage of the 
light may occur. We do not deal with this explicitly, but we note that this can generally be 
achieved by choosing the reduced optical path length large enough. (See Section 2 for the 
definition of the reduced optical path length.) 



2 



the unit sphere S'^, and T is contained in a plane perpendicular to the direction 
of the output beam. (See again Figure 1.) The input and output intensities 
are given by nonnegative functions I{m),m G D, and L{x), x eT. 

We represent the two reflectors via the polar radius p{m),m e D for the 
flrst one, and as the graph of a function z{x), x e T for the second one. See 
again Figure 1. One of the main results is that finding the functions p{m) and 
z{x) is equivalent to solving the following constrained optimization problem 
for certain transforms of p(m) and z{x): 

Minimize / log p{m) I (m) da + / log z{x)L{x)dx (1) 



subject to logp(m) + logz{x) > log K{m,x) for all m E D,x E T. (2) 

Here logp(m) and logz{x) are certain transforms of p{m) and z{x), respec- 
tively, given explicitly in Definitions 2.1 and 2.2 below. The constraint is given 



Here m — {mx,mz) eDCS'^,xefc R^, and £ > is a constant. For 
details see the text below. 

Note that this reformulation not only gives a theoretical existence result for 
solutions, but it also translates into a practical method for finding the solution. 
In fact, the discretization of the constrained optimization problem (l)-(2) is a 
standard linear programming problem and can be solved numerically. 

As mentioned before, the approach of this paper has a strong connection 
to, and is in fact motivated by, the theory of optimal transportation. (See 
for example [1, 2, 3], and in general the recent survey [10] and the extremely 
extensive bibliography cited there.) Consider the ray tracing map, or reflector 
map, 7 : D — > T. So a ray emitted in the direction m E D will be transformed 
by the reflector system into a ray labeled hy x — 7(m) E T. (See again 
Figure 1, and also Figure 4.) 

Consider the optimal transportation problem for moving the measure I{m)da 
on D Q to the measure L[x)dx on T C via a transformation P : D ^ T 
in such a way that the "transportation cost" 



is maximized'^. We prove that the ray tracing map 7 actually solves this 
^It is unclear to the author whether the cost C(P) has any physical meaning. If max- 



by 



log i^(m, x) = log 



£-{ma;,x) 1 



2£(£2 



x\^){l + m,) 4^2 
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problem. This gives a uniqueness result for the reflector map, see Theorem 
7.3. 

The structure of the paper is as follows: In Section 2, we state the problem 
and our notation precisely and introduce the transforms p{m) and z{m). In the 
next section, we give certain analytical definitions of our concepts of reflectors 
and the corresponding reflector map. In Section 4, we then show that these 
analytical constructions correspond to some geometric constructions. This 
geometric content is that we seek to represent the first reflector as the boundary 
of the intersection of a certain family of spheroids, and we seek to represent the 
second reflector as the boundary of the intersection of a family of paraboloids. 
We also justify that our abstract deflnition of the reflector map 7 is consistent 
with the optical definition. In Section 5, we formulate the rcfiector problem 
and show that it is equivalent to the constrained optimization problem (l)-(2). 
Then in the next section, we prove the existence of solutions, which follows 
from a standard compactness argument. Finally, in Section 7, we show the 
connection to the transportation problem mentioned above, and finally state 
the main theorem, Theorem 7.3, on the existence and uniqueness of solutions 
to the reflector problem. 

Throughout, the treatment is very similar to the papers [5,4], which treated 
similar but distinct beam shaping problems. The main difference is that the 
problem at hand requires the introduction of the transforms p(m) and z(x), 
and a more complicated form for the cost function log K{m,x). These com- 
plications arise mostly because the flrst reflector is best described using polar 
coordinates, and the second reflector is best described using Cartesian coordi- 
nates. For the sake of completeness and being self-contained, we include most 
of the proofs here, leaving out some details if they can easily be filled in from 
[5] and [4]. 

2 Statement of the problem and assumptions 

We first fix our notations and assumptions in this section. Consider the con- 
figuration show in Figure 1. A point source located at the origin O — (0, 0, 0) 
generates a spherical wave front over a given input aperture D contained in 
the unit sphere S"^. The latter is required to be transformed into a beam of 
parallel rays propagating in the direction of the negative 2;— axis by means of 
a system of two reflectors. A cross section of the output beam is specifled as 

imizing a transportation cost seems unintuitive, an equivalent formulation is of course to 
minimize —C{P) 
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a set on a plane perpendicular to the direction of propagation. Such a cross 
section is called the output aperture, and denoted by T. 

We denote points in space by pairs {x, z), where a; G is the position 
vector in a plane perpendicular to the direction of propagation and 2; e R is 
the coordinate in the (negative) direction of propagation. See again Figure 1 
for our convention on the direction of the 2;— axis. Points on the unit sphere 

will typically be denoted hj m E S"^; their components are also written as 

= i'tTT'x, ^z) with |ma;p + ml = 1. 

We fix the output aperture in the plane z — —d. We will seek to represent 
the two reflectors as the graph of the polar radius p{m) and as the graph of a 
function z{x), respectively, as shown in Figure 1. That is 

Reflector 1: = {pim) ■ m \ m E D}, 
Reflector 2: F^ = {{x, z{x)) \ x ef}. 

We now have the following assumptions on the input data. Let D be 
an open, nonempty subset of with closure D such that (0, 0,1) ^ D and 
(0, 0, —1) ^ D. So there is some £ > such that for m — {rrix, niz) G -D, we 
have 1 — e > rriz > — 1 + e. Let further T to be an open, bounded, nonempty 
subset of M^. Denote its closure by T. 

The geometrical optics approximation is assumed. It follows from general 
principles of geometric optics that all rays will have equal length from (0, 0, 0) 
to the plane z — —d; this length is called the optical path length and will be 
denoted by L. We define the reduced optical path length as i — L — d. In terms 
of the mathematical problem, £ is an input parameter. We pick the reduced 
optical path length i > large enough so that the following conditions are 
satisfied: 

i > max \x\, and (3) 
^ ^ - — >OiorallmeD,xef. (4) 



(£2- |x|2)(l + m^) 2£ 

(It is not hard to see that the second condition is indeed satisfied for large 
enough £.) 

We'll also use the convenient notation 

Finally, we define the following two transformations, which are central to the 
analysis: 
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Definition 2.1. Let z = z{x) he & continuous function defined on T C R^. 
Then define the function 

z{^)-S--^^^ iorxef. (5) 

Definition 2.2. Let p — p{m) be a continuous function defined on D C S'^ 
with p > 0. Then define the function 

p(m) = -S + — — r-^ — - for m e D. (6) 

These transformations are obviously invertible, namely 

z(x) = {-z{x) + 5){f - \x\^) for xef, 
p(m) ^ {2(m^ + l){p{m) + 5))'^ formed. 

The following lemma is now obvious with the above formulas. 

Lemma 2.3. (i) The transformation C{T) C{T),z ^ z, is a bijection. 

(a) Denote by C>q{D) the set of all positive continuous functions on D. 
Then the transformation C^q{D) — > {ct e C{D) | a > — ^ p, is a 
bijection. 



3 Reflector pairs and the reflector map 

We now come to the central definition of this paper, namely that of a reflector 
pair. We first give the analytic definition. In the next section, we look at the 
geometric interpretation. In preparation, we first define the following function: 

Definition 3.1. Define the function 

K{m, x) = S ' J:' - 5^ for m = (m^, m^) e D,x eT. 

[i^ — + rriz) 

Note that K{m, x) > by (4). The function log7^(m, x) can be interpreted as 
some kind of cost of transporting a unit of energy from m to x. (See Section 7 
for more details.) 

We are now ready to define the notion of a refiector pair. 
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Definition 3.2. A pair 
p,z > and 

p{m) - 
z{x) - 

Here we used the definitions of p and z from Definitions 2.2 and 2.1. Note 
that the suprema on the right hand sides are in fact attained. Also note that 
for a reflector pair (p, z) , we have 

p{m) z{x) > K{m, x) for allm e D,x e f. (9) 

We will show later that reflector pairs are in fact concave (see Propo- 
sition 4.8) and that logp and logi are uniformly Lipschitz continuous (see 
Proposition 6.2). 

Finally, we define the reflector map, or ray-tracing map, associated to a 
reflector pair. Again, the choice of terminology will become clear when we 
consider the problem from a geometric viewpoint in Section 4.4. 

Definition 3.3. Let {p,z) e C>o(-D) x C{T) be a reflector pair. Deflne its 
reflector map, or ray tracing map, as a set- valued map 7 : — > {subsets of T} 
via ^ 

7(m) = {x e T I p(m) = K(m,x)} for m E D. 

z{x) 

Clearly 7(771) 7^ for any m E D hy (7). We will show later that 7(m) is in 
fact single-valued for almost all m E D. (See Proposition 4.10.) We may thus 
regard 7 as a transformation 7 : D ^ T. 

4 Reflector pairs: Geometric viewpoint 

In this section, we will investigate the definition of refiector pairs from a ge- 
ometric point of view. We show that the reflectors can be obtained as the 
boundary of certain convex sets. These sets in turn are the intersections of a 
family of spheroids (for reflector 1) and paraboloids (for reflector 2), respec- 
tively. 

In the first two sections, we will use the following notation, the pointwise 
analogue of Definitions 2.1 and 2.2: For x e and z > 0, and m E and 



(p, z) e C>o(-D) X C{T) is called a reflector pair if 



sup ( ■^^K{m,x) ) for m G D, (7) 



z{xy 



sup ( -z^K(m,x) ] for a; e T. (8) 
men \p{m) 
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Figure 2: LEFT: The set Sx,z is a spheroid with foci at the origin O and at 
(x, z). Sx,z is the convex set bounded by S^^z- The sketch shows a cross section 
through a plane containing the axis defined by the foci. The sets S^^z and S,j:^z 
are obtained by rotating the gray elhpsc around the axis through the foci as 
indicated. RIGHT: The set Px,z is a paraboloid with focus at p • m and axis 
parallel to the the convex set bounded by Pp.m- The sketch 

shows a cross section through a plane containing the axis and the focus. The 
sets Pp.m and Pp.m are obtained by rotating the gray parabola around the axis 
as indicated. 

p > 0, write 

4.1 The spheroids Sx,z 

We define first a family of spheroids, indexed by points (x, z). 
Definition 4.1. Let x ef and 2; e M such that z > 0. Define the set 

Sx,z^ {p-m\me S'^,p^]:K{m,x)}. (12) 

z 

Here we used the notation from (10) and (11). 



Lemma 4.2. Using the notation from the previous definition, the condition 
z > implies i — z > \{x,z)\. (Here |.| denotes the standard Euclidean vector 
norm.) The set Sx,z is given by the following equation for p > 0, m & : 

p + \{x, z) — p ■ m\ = i — z. (13) 

Geometrically, Sx,z is a spheroid whose foci are the origin O — (0, 0, 0) and 
the point (x, z). (See Figure 2.) 

Proof. Note first that the condition 5 > imphes 



z = {5- ~z){f - \x\^) < 5{t - \x\^) = -- 
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and so 




> 0. 



Also, 



It follows that 



i- z> \{x,z)\. 



A straightforward but lengthy algebraic computation now yields that the con- 
dition in (12) is equivalent to 



Here p = i^im^ + l)(p + b))'^ > 0, as p > 0. We have 

I — z > \{x^ z^\'>i p — \ {x^ z^ — p ■ m\ — p — \i — z — 

and thus £ — z — p > — 1£ — z — p\. It follows that 

i — z — p — \i — z — p\. 

This together with (14) gives equation (13). Now consider the geometric con- 
tent of equation (13). We can immediately read off that the set S^^z consists 
of all points that satisfy that the sum of the distances to the points (0, 0, 0) 
and (x, ^) equals i — z. (See Figure 2. Also note that i — z > \{x,z)\ implies 
that the set Sx,z is nonempty.) This is by definition a spheroid. □ 

Definition 4.3. Let as before a; G and z G M such that z > 0. Denote by 
^^,2 the closed convex set bounded by Sx,z- Thus Sx,z is given by 



4.2 The paraboloids Pp.^ 

Similarly to the previous section, we now define a family of paraboloids Pp.m 
Definition 4.4. Let m E D C S'^ and p > such that p > 0. Define the set 



{£ — z — pY = \{x, z) — p ■ m 



(14) 



S{x,z) = {p • m I m e 5*^, p > -K{rn, x)}. 



(15) 



p-m 



{{x,z) I x e 



-K{m,x)}. 



(16) 



P 



Here we used again the notation (10), (11). 
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Lemma 4.5. With the notation of the previous definition, consider new coor- 
dinates {p, q) obtained by a shift by p ■ m, that is (p, q) = {x, z) — p ■ m. Then 
the set Pp.m is given by the equation 



4ag = IpI^ - 4a^ (17) 

where 

A5 

Geometrically, Pp.m is a paraboloid with focus at p ■ m and axis parallel to the 
z—axis. It opens in the directions of the negative z—axis. (See Figure 2.) 

Proof. A straightforward computation yields that (16) is equivalent to (17). 
Note that a paraboloid with focus at (0, 0, 0) and focal parameter 2a has the 
equation Aaz = \x\'^ — Aa^. Also note that a < follows from p > 0. □ 

Definition 4.6. Let as before m E D C S"^ and p > such that p > 0. Denote 
by Pp.m the closed convex set bounded by Pp.m- Thus Pp.m is given by 

Ppm = {{x,z) \x eR^,z> \K{m, x)}. (18) 



4.3 Geometry of reflector pairs 

We now investigate the geometric content of the definition of reflector pairs. 
For this, suppose that (p, z) e C>o(i)) x C{T) is a reflector pair. Consider the 
two sets 

S^f]S,,,^,), (19) 

xef 

mG-D 

Note that none of the sets Sx,z{x) and Pp(m)m ^-re empty as in fact p(m) > 
for all m G -D by definition of refiector pairs and z{x) > for all a; G T. A 
compactness argument yields that P and S are nonempty as well, and they 
are convex sets, since they are the intersections of convex sets. Moreover, the 
algebraic representations (15) and (18) along with the definition of refiector 
pairs immediately give rise to the following geometric facts: 

Proposition 4.7. Let {p, z) G C>o(-D) x C{T) be a reflector pair. 
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Figure 3: Geometry of reflector pairs. LEFT: The graph Tp of p(m) is obtained 
by intersecting the boundary of the convex set S with the cone M+D. RIGHT: 
The graph of z{x) is obtained by intersecting the boundary of the convex 
set P with the cyhnder T x R. 

(i) The graph of p, 

Fp = {p{m) ■ m \ m & D}, 

is contained in the boundary dS of the convex set S given in (19). In 
fact, 

where ■ D — {r-m | r>0,me D} is a cone with vertex at the origin 
and cross section D. (See Figure 3.) 

(a) The graph of p, 

V, = {{x,z{x)) \xeT} 

is contained in the boundary dP of the convex set P given in (20). In 
fact, 

F^ = aPn(TxR), 

where T x R = {(x, /i) | x e T, /i e R} is the cylinder with base T. (See 
Figure 3.) 

We may thus think of the two reflectors as the "envelopes" of certain fami- 
hes of spheroids and paraboloids, respectively. An immediate corollary is also 
the following result: 

Proposition 4.8. // {p,z) G C>o(-D) x C(T) is a reflector pair, then z is a 
concave function and p is the radial function of a convex set. In particular, both 
p and z are locally Lipschitz continuous and almost everywhere differentiate. 
(Here "almost everywhere" refers to the standard Lebesgue measures on 
and R^, respectively.) 

Proof. The first part of the above statement follows immediately from the pre- 
vious discussion. The Lipschitz and differentiablility properties are standard 
results from convexity theories, see [9], Theorem 1.5.1 and Section 1.7. □ 
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4.4 Geometry of the reflector map 



We now investigate the geometry of the reflector map. For this, the following 
terminology is useful: Suppose (p, z) e C>o(-D) x C{T) is a reflector pair. For 
m G -D, if p(m) ■ m G Fp fl S^^zix), say that the spheroid Sx,z{x) is supporting 
to the graph Tp at the point p(m) ■ m. 

Similarly, for x G T, if {x, z{x j) G F^ n Pp{mym, say the paraboloid Pp{mym 
is supporting to F^ at (x, z{x)). 

Note that Sx,z(x) is supporting to Fp at the point p{rn) ■ m if and only if 
Pp{m)-m is supporting to F^ at {x,z{x)). This is because both statements are 
equivalent to 

z{x) p{m) = K{m, x) 

by the deflnitions (16) and (12). 

We also have the following geometric lemma: 

Lemma 4.9. Let Xi,X2 G T be two distinct point: Xi X2- Suppose the 
spheroids Sx^^z{xi) ond Sx2,z(x2) ^'^^ supporting to Tp at the same point p{m)-m. 
Then Sx^x^i) and Sx2,z{x2) intersect tranversally at p{m) • m. 

Proof. Assume the contrary, that is, that Sx^^z{xi) and Sx2,z{x2) intersect tan- 
gentially to each other. Since the two spheroids share the focus O, it follows 
from basic properties of ellipsoids that the two line segments (xi, z{xi)), p{m) ■ m 
and (a;2, ^(2^2)), p{fn) • are parallel. Thus the three points (xi, z{xi)),{x2, -2(2:2)) 
and p(m) ■ m are collinear. On the axis through these three points, the point 
p{m) -m cannot lie between the other two points; that is p{m)-m must be one of 
the end points of the hue segment defined by the three points. But this contra- 
dicts that the paraboloid Pp[m)-m contains both (xi, z{x\)) and (^2, z{x2)). □ 

Now consider the reflector map 7 associated with a reflector pair (p, z) 
as defined in Definition 3.3. In the language defined above, we may now 
say that 7(m) is the set of all points x E T such that that the spheroid 
Sx,z{x) is supporting to the graph Fp at the point p{m) ■ m. Consider the 
case where for some m E D, the set j{m) contains more than one point, say 
{xi,X2} C 7(m). By the previous lemma, the two spheroids S^j^^zixi) and 
Sx2,z{x2) intersect transversally at p{m) • m. Thus p is not differentiable at 
m. But the set of such points has measure zero by Proposition 4.8. We thus 
immediately have the following result: 

Proposition 4.10. Let {p,z) G C>o{D) x C{T) be a reflector pair. Then its 
reflector map 7: D ^ T is almost everywhere single-valued. That is, the set of 
points m where 7(771) is not single valued has measure zero with respect to the 
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Figure 4: Geometry of the reflector map: A ray emitted from the origin O will 
be reflected to a ray traveling in the negative direction labeled by x = ^{x). 
See the text for details. 

standard measure on as a submanifold of the measure space (M^,/x), where 
II is the standard Lebesgue measure. 

Remcirk 4.11. We now justify the terminology of the reflector map from 
an optical point of view. See Figure 4 for the following considerations. Let 
{p,z) G C>o(-D) X C{T) be a reflector pair. Let m E D such that 7(m) is 
single valued. We show that under the geometric optics approximation, a ray 
emitted in the direction m E D will be reflected off the first reflector Fp and 
then the second reflector F^ in such a way that the reflected ray is parallel to 
the negative 2;— axis and that it intersects a plane perpendicular to the 2;— axis 
in the point x = 7(772). 

Consider the reflection off reflector 1 first. Since the spheroid Sx,z{x) is 
tangential to the reflector Fp, the ray will be reflected off Fp the same way it 
would be reflected off S^Xx)- By the geometrical properties of spheroids, this 
means that the ray is reflected towards the focus {x, z{x)). There, the ray will 
encounter F^. It will be reflected the same way as it would be reflected off the 
paraboloid Pp(m)-m, that is, in the direction of the negative 2;— axis. 

Thus our definition of the reflector map is in agreement with the physical 
law of reflection. In the case when 7(m) is multi-valued, the first reflector has 
a singular point and a ray will split up into a cone of light rays. These rays 
will generate a set of directions whose projection onto a plane perpendicular 
to the 2;— axis is 7(m). This is consistent with the physical phenomenon of 
diffraction at singularities. 

The following statements about the reflector map 7(m) are analogous to 
Theorem 4.8 and Lemma 4.9 in [5]. See this paper and the the reference [7] 
for further details on the proofs. 

Theorem 4.12. Let B denote the a algebra of Borel sets on T. Then for any 
subset T & B, 7~^(t) is measurable relative to the standard Lebesgue measure 
of D. In addition, for any non-negative locally integrable function I on D, the 
function 
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is a non-negative completely additive measure on B. (Here da is the standard 
measure on S"^ C M^.j 

Lemma 4.13. With the notation of the above Theorem, let h be a continuous 
function on T. Then we have 



5 The reflector problem and an equivalent con- 
strained minimization problem 

Let now / be a nonnegative, integrable function on D, and L be a nonnegative, 
integrable function on T, such that 



We may interpret / and L as the intensity distribution functions of the hght 
beams on the input and output apertures, respectively. The above integral 
condition is simply (total) energy conservation. 

In this section of the paper, we now formulate the reflector problem. More 
specifically, we can call this formulation a "weak" version of the reflector prob- 
lem since we do not require the input functions / and L to be differentiable, 
nor do we require the reflectors to be smooth surfaces. 

We then formulate a second problem, which is an inflnite dimensional linear 
programming problem. One of the main results is that the two problems are 
in fact equivalent. This is stated and proved at the end of this section. 

Let us flrst formulate the reflector problem: 

Problem I. (Reflector Problem) For given input and output intensities I 
and L satisfying (21), find a pair {z,p) e C>o{D) x C(T) that satisfies the 
following conditions: 

(i) {z,p) is a reflector pair 

(a) The reflector map 7: D — > T satisfies 





(21) 




for any Borel set t CT. 
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This formulation builds on the geometrical interpretation of reflector maps 
as presented in Section 4. Note that condition (ii) is local energy conservation. 
We have the following immediate corollary from Lemma 4.13: 

Corollciry 5.1. Let {z,p) be a solution to Problem I. Then we have 

/ h{'-){m))I{m)da = I h{x)L{x)dx 

JD JT 

for all functions h e C{T). 

Before we now formulate Problem II, we define the following function space: 
Definition 5.2. Define the set of admissible functions as 

Adm(L>,r) = {(r,C) e C{D) x C{f) \ r{m) + Ci^) > log K(m,x) 

for aWm e D,x e f}. 

Problem II. Minimize the functional 

T{r, C) = / r{m)I{m)da + / Cix)L{x)dx 

JD JT 

on the space Adm{D,T). 

The two problems are equivalent, as expressed in the following theorem. 

Theorem 5.3. Let {p,z) G C>o(-D)xC(T) be a reflector pair. Then (log p, log 5) e 
Adm(Z),T'). The following statements are equivalent: 

(i) {p,z) solves the Reflector Problem L 

(ii) (logp, logS) minimizes the functional on Adm{D,T). 

Proof. The statement (log p, log z) e Adm(L), f) follows immediately from (9). 
The proof of the equivalence of (i) and (ii) is analogous to those of Theorem 
5.2 in [5] and Theorem 3.4 in [4]. See also Theorem 1 in [3]. Since this theorem 
is central to this paper, we give an outline of the proof, omitting some of the 
technicalities, which can be filled in with the above references. 

(i) ^ (ii) : Suppose (p, z) solves Problem I with corresponding reflector 
map 7. Let (r, ^) e Adm(D,T). Then for any m & D such that 7(m) is 
single-valued, we have 

r(m) + C(7(^)) ^ logi^'(m, j{m)) — log p(m) + log5(7(m)). 
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This yields 



r{m)I{m)da+ / ({■j{m))I{m)da 



D 



> 



I log p{m) I {m) da + / log z {'y{m)) I {m) da. 

D J D 



Now using Corollary 5.1 gives ^(r, C) > ^(logz, logp). 

(ii) =^ {i) : The main idea is that the Euler-Lagrange equations of mini- 
mizing !F are equivalent to the equality 



for all functions h G C{T). This implies that (p, z) solves Problem 1. 

Let thus h G C{T). Let £ > be a small parameter. To bring out the main 
ideas, we present a formal calculation, assuming expansions in e are valid. A 
completely rigorous treatment is possible; indeed the proof in [5] can easily be 
modified to the problem at hand. 

Define perturbations {reXe) ^ Adm{D,T) of (log p, log 5) via 



Let now Xe be a point where the supremum in the definition of rs{m) is at- 
tained. Expanding x^ in e yields 



Thus again an expansion in e gives 

r^{m) = -Cs{xe) + log K{m,Xe) 

= — log z(7(m)) + log K{m, 7(m)) — e h{'^{m)) -\- 0{e^) 
= logp(m) - eh{^{m)) -\- 0{e^). 

Thus, using the fact that (log p, log z) minimizes 




Ce{'^) — log z{x) -\- e h{x) for x e T 

re{m) — sup {—Q{x) + log Kim, x)) for m & D. 
xet 



Xe^-f{m)+0{e). 




de £=0 



j^{re,Ce)^- I h{j{m))I{m)da+ / h{x)L{x)d: 



Ix. 




This completes the sketch of the proof. 



□ 
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6 Existence of solutions 

Note that Theorem 5.3, while showing that Problems I and II are equivalent, 
does not state that solutions exist. We prove this in the present section. 
In the following, fix some point m* & D. 

Definition 6.1. Set 

Admrefl(5,f ) = {(logp,log5) e C{D) X C{f) I 

(p, z) is a reflector pair with logp(m*) = 0}. 

Note that Adm,efl(£', f) C Adm(L>, f). 

Proposition 6.2. The family of pairs of functions Admrefi(5, is uniformly 
Lipschitz continuous in each entry. That is, there are constants Ki,K2 > 
such that 

\C{xi) - C{X2)\ < Ki ■ \xi - X2\ 

\r{mi) —r{m2)\ < K2 ■ d{mi,m2) 

for all {r,() G Admrefl(-D, T), xi,X2 G T, mi,m2 G D. Here d{mi,m2) = 
dist52 (mi, 777.2) is the intrinsic distance on S^. The constants Ki,K2 only 
depend on D, T; explicitly, we have 

Ki — rnax _ | V^; log K{m, x) \ , 

meD,xeT 

K2 = rnax _ | log K{m,x)\. 

meD,xGT 

(Here Vm denotes the gradient with respect to the variable m on the sphere 
S'^, and Vx is the gradient with respect to x.) 

Proof. We prove the inequality for (; the proof of the other inequality is com- 
pletely analogous. Let Xi,X2 G T and {r,() G Admrcfl(-D, T). Assume ({X2) < 
C{xi), otherwise relabel Xi,X2. Let mi G -D be such that Ci^i) + ^"(^i) = 
\ogK{mi,Xi). (Such an mi exists by (8).) Then 

< Ci^i) - ({^2) = -({^2) - rimi) + log K{mi, Xi) 
< logK{mi, Xi) — log K {mi, X2) 

<( max _\VxlogK(m,x)\] ■ \xi — X2\. 

\m£D,xeT J 

□ 
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Proposition 6.3. Functions m Admrefl(I?, T") are uniformly bounded. Specif- 
ically, we have 

\r{m)\ < Ki diamp) 

\C(x)\ < max|logX(m*,x')| +i^2 diam(r). 

x'eT 

for all (r, () G Admrefl(-D, T), x & T, m & D . Here Ki, K2 are as in Proposi- 
tion 6.2, and diam(T) = max^^ •j.^gf a;2| an(i diam(Z)) — maXj^^ ^j^g^ d(mi, 7712) 
are the diameters ofT and D, respectively. 

Proof. Let (r, C) e Admrefl(-D, T), x E T, m E D. By Proposition 6.2, we have 

|r(m)| < |r(m*)| + Kid{m,m*) < Ki diam(Z)). 

Also, for some x* G 7(771*), we have C,{^*) — log i^(m*, x*), and thus again by 
Proposition 6.2, 

\C{x)\ < \C(x*)\+K2\x-x*\ < max|logir(m*,x')| +i^2 diam(f). 

x'eT 

□ 

We can now prove the existence of a solution to Problem II. With Theo- 
rem 5.3, this immediately implies the existence of a solution to the reflector 
problem. 

Theorem 6.4. The functional attains its minimum on Adm{D,T) . More- 
over, this minimum is actually attained at some {r,Q G Admi.efl(5, T). 

Corollary 6.5. The Reflector Problem I has a solution. 

Proof. By Propositions 6.2 and 6.3 and the Arzela-Ascoli Theorem, the func- 
tional T attains its minimum on Admrefl(-D, T). 

We now show that this minimum is also the minimum of J- on the larger 
set Adm(L>,f). For this, let (r, C) e Adm(L>,f ). Define 

r*(m) = sup{— + log K{m, x)} ioi m E D 
xef 

C*{x) — sup{— r*(77i) -\-log K{m,x)} for x G T 

mG-D 

It follows that 

r*(m) < r(m), C*(^) ^ C(^) 
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for all m e D, X e T. Note that C*(a;) < C,{x) implies 

sup{— C*(a:;) + \ogK{m, x)} > sup{—C{x) + log K{m, x)} — r*{m). 

But —C*{x) + log K{m, x) < r*{m) for all x e T, m e D implies 

sup{— + log K{m,x)} < r*{m). 

xef 

It follows that 

r*(m) = sup{— + log K{m,x)} for m & D. 

Thus {r*X*) = 0-Og p,logz) for some reflector pair {p,z). Note that 

{r* - r*{m*), C + r*{m*)) e Adm,^Q{D, f), 

where the left hand pair denotes functions shifted by the constants ±r*(m*). 
Using the above, we now have 

J^{r* -r*{m*),C + r*{m*)) = J^{r*,C) = [ r*{m)I{m)da + j_C{x)L{x)dx 

JD JT 

< j r{m)I{m)da+ / C{x)L{x)dx ^ J^{rX)- 

JD JT 

This shows that indeed the minimum of T on Adm(Z),T) is attained on 
Admrefl(-D, T). The corollary is an immediate consequence of Theorem 5.3. □ 



7 A uniqueness result and connection to an 
optimal transportation problem 

In this section, we show the connection to an optimal transportation problem. 
This is again quite analogous to the reflector design problems in [5] and [4]. 
This connection allows us to formulate a uniqueness result for the Reflector 
Problem I. 

For the formulation of the problem, we need the concept of a plan in this 
context: 

Definition 7.1. A plan is a map P: D ^ T that is measure preserving, that 
is, we have 

/ h{P{m))I{m)da = I h{x)L[x)dx 

J D JT 

for any function h e C{T). 
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A plan P needs only be defined almost everywhere on D. 

The optimal transportation problem associated with intensities /(m) and 
L{x) satisfying energy conservation (21) and with cost function log i^(m, a;) is 
the following: 

Problem III. Maximize the transportation cost 



among all plans P. 

This problem is again solved by the reflector map of a solution to the 
Reflector Problem I. In fact, we have the following theorem: 

Theorem 7.2. Let {z,p) be a solution to the Reflector Problem I. Then the 
corresponding reflector map j is a plan, and it maximizes the transportation 
cost (22) among all plans. Any other cost maximizing plan is equal to 7 almost 
everywhere on supp / \ {m e D | I(m) — 0}. 

Proof. Let {z, p) be a solution to the Reflector Problem I with corresponding 
reflector map 7. By Corollary 5.1, 7 is in fact a plan. Now let P be another 
plan. Then 

logp(m) + log 5(P(m)) > \ogK{m,P{m)) 
for almost all m E D, and equahty holds iff P(m) = 7(m). Thus 



Thus 7 indeed maximizes the transportation cost among all plans. Moreover, 
if equality holds in the above estimate, then 7(771) = P(m) or /(m) = for 



P^C{P)= / log K {m,P{m)) I {m)da 



(22) 



D 




almost all m e D. 



□ 
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We can now state the following theorem, summarizing the main result 
on the existence of solutions to the Reflector Problem I with an additional 
uniqueness result: 

Theorem 7.3. (Existence and uniqueness for solutions to the Reflector Prob- 
lem I) There exist solutions to the Reflector Problem I. If (p, z) is a solution, 
then both p and z are Lipschitz continuous. The corresponding reflector map 
^y: D ^ T is single valued almost everywhere on D. If {p,z) and {p',z') are 
two solutions with reflector maps 7 and 7', respectively, then 7(m) = 7'(?7t.) 
for almost all m e supp / \ {m e D | /(m) = 0}. 

Proof. The existence of solutions to Problem I was already obtained in Corol- 
lary 6.5. If now 7 and 7' are two reflector maps corresponding to two solutions, 
then both maximize the transportation cost C by Theorem 7.2, and hence we 
have 7(x) = 7(0;') for almost all supp / \ {m e D | /(m) = 0}. □ 

There is a number of open questions for further investigations. For in- 
stance, we have assumed certain constraints on the aperture D, in particular 
(0,0,-1) ^ D. From the physical intuition about the problem, these con- 
straints appear to be unnecessary. It would be interesting to extend the theory 
to these cases as well. Furthermore, a further exploration of the regularity of 
solutions (p, z) depending on the intensities / and L would be desirable. It is 
expected that some of the currently rapidly growing research on optimal trans- 
portation may carry over here. (See the recent survey [10] and the extremely 
extensive bibliography cited there.) 

Acknowledgement: The author would like to thank V. Ohker for useful 
discussions. 
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